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Abstract 



We evaluate the contribution of the final state interaction (FSI) due 
to single pion exchange inelastic scattering for D + — > K°*tt + and 
D + — > K°p + processes. The effects are found to be very significant. 
The hadronic matrix elements of the weak transition are calculated in 
terms of the heavy quark effective theory (HQET), so are less model- 
dependent and more reliable. 



PACS number(s): 13.75.Lb, 13.20. Gd, 13.25, 14.40 



I. Introduction 

To understand precisely the mechanism governing weak transition processes 
where fundamental physics and possible new physics apply, one needs to face a 
synthesis problem including evaluation of the hadronic matrix elements and final 
state interactions (FSI). Unless one can fully understand the side effects, such as 
FSI, one can hardly extract information about new physics correctly. At least one 
needs to know the order of magnitude of FSI and determine its significance. 

Due to the success of the Standard Model, the Hamiltonian for the weak transi- 
tion c ^ s + u + d is well understood |]T], and thanks to the heavy quark effective 
theory (HQET) ||, we have some reliable ways to handle the hadronic matrix ele- 
ments for the b — > c transition. The developments of the chiral lagrangian || enable 
us to estimate matrix elements from heavy mesons to light pseudoscalars. According 
to the recent work by Roberts and Ledroit [[|, the transition matrix elements from 
B, D to K(*\ become calculable in a unique theoretical framework. 

On the other hand, to really test the theory concerning matrix elements and 
HQET, one needs to calculate FSI effects which sometimes become very significant 
or sometimes can be negligible. 

In the early work, Buccella et al. found that the calculated branching ratios 
for D + — > K*°7T + and D + — > K°p + deviated significantly from the data and even 
using recent data, the difference is still obvious. We use the HQET to recalculate 
(see below for details) and find that the discrepancy with data is reduced, but still 
exists. This disagreement may be caused by the final state interaction (FSI). In fact, 
many authors have studied the problem of FSI in some D and B decays [|5], [| [7], § . 
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The FSI due to s-channel resonances were found to be very important for some final 
channels |5|, ||. However for the processes D + — > K* tt + and D + — > K°p + , the final 
states have isospin 1=3/2, and therefore have no s-channel resonance FSI. Zheng @ 
calculated the elastic FSI effects in B — > DK caused by t-channel meson exchange, 
and obtained a very small phase, so it does not make a substantial effect to the 
measured rate. Very recently, Donoghue et al. indicate that inelastic scattering 
may dominate the FSI ||. We concur with this. Our previous work in \l/ decays and 
pp annihilations indicate that when a single— 7r exchange in the t-channel can be 
realized, the corresponding mechanism would make a significant contribution to the 
FSI and would very probably dominate. 

For D + — > K*°7i + and D + — > K°p + , the dominant FSI should be the inelastic 
rescattering between K*°tt + and K°p + , as shown in Fig.l. In this work we estimate 
this inelastic FSI effect. The hadronic matrix elements are evaluated in terms of the 
method given in ref.PJ which can alleviate the model-dependence of the calculations. 

In Sec. II, we give the formulation of the transition amplitude with and without 
considering the inelastic scattering of K*°7t + <-> K°p + , and in Sec. Ill, we present our 
numerical results, while the last section is devoted to our discussion and conclusion. 



II. The formulation 

In the weak interaction, isospin is not conserved. There are four possible VP 
decay modes for D + , i.e., K*°n + , Kp + , K* + n° and K + p°; among them I 3 (K*°n + ) = 
h(K°p + ) = 3/2 and I 3 (K* + p°) = I 3 (K* + tt°) = 1/2. However, from the quark 
diagrams |10| , one finds D + — > K* + tt° or D + — > K + p° can only be realized via a 



Cabibbo doubly suppressed channel, so must be very small and can be neglected. 



The FSI conserves isospin. The K*°tc + and K°p + cannot rescatter into K* + p° and 
K* + 7i°; while K*°tt + <-> K°p + can be realized by the t-channel single pion exchange 
diagrams shown in Fig.l and can be very important. Therefore, we only need to 
consider production of D + —>■ K*°tc + and D + — > K°p + , as well as their mutual 
conversion through inelastic scattering, 
(i) Without the FSI. 

The weak interaction hamiltonian for non-leptonic decay is given as 
G 

H w = -Lv ud V^[ Cl s ltl {l - 75)^(1 - ls)d + c 2 s 7m (1 - 75)^(1 - 75 )c] (1) 

where V u d and V cs are the Cabibbo-Kobayashi-Maskawa entries. The reactions con- 
cerned are Cabibbo favored processes. The color indices are dropped out as well 
understood, and the coefficients c\ and c 2 are obtained from the renormalization 
group equation [IJ. 

In the calculations, we use vacuum saturation and ignore the W-exchange (an- 
nihilation) contributions JTU|]. The transition matrix elements for D + — > K*°tt + and 



D + -> #V are : 



Tl (D + -> i?*°7r+) = ( Cl + + 5)c 2 ) < tt + |«7 m (1 - 7s)rf|0 >< K*°\s llA (l - l5 )c\D + > + 

+ (ca + ^r(l + <J) Cl ) < iT |s 7 M (l - 75)^10 >< tt + |^7 m (1 - T5 )c| J D + > 

= -(ci + ^(1 + 5)c 2 )/^Tr{[(£ 3 + ^ 4 )e A „ ■ v 

+ £. + Ac&)mi - 75) ■ - m 



+ (ca + ^(l + *)d)/*.C™x. I V{(ei + ^6)757^ ■ ^^75(1 - j)}, 



(2) 
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and 

T 2 (D + ^K°p + ) = ( Cl + ^(l + 5)c 2 )<p + \uY(l-l 5 )d\0><K°\s lfl (l- l5 )c\D + > + 

+ (c 2 + ^-(1 + S) Cl ) < K°\s^{l - 75 )d|0 >< p + \u ltl {l - ^)c\D + > 

= (ci + ^(1 + S)c 2 )f p m p e^Tr{(^+^ 2 ) l5l , ■ ^^75(1 - fl} 

- (c 2 + i(l + S) Cl )f K p^ K Tr{[(^ +^e*)e p • « 



+ ^5+M6)]7m(1-75)-2\/ M ^(1-^}, (3) 

where u is the four-velocity of Z} + as = M^v^ and & (i=l,...,6) are functions 
of momenta given in ref.Q; 5 is a non-factorization term and cannot be evalu- 
ated in perturbative QCD [[11], |I2 . Recently, Sharma et al.p^| investigated the 



non- factorization effects in D — > PV^ decays. Blok and Shifman||14|| give a more 
theoretical estimation of the factor as 



xml H 



where x ~ 1 and mo-// is a numerical factor. Admittedly ||1 4j| , one cannot take the 



number as accurate, so we keep it as a parameter in the region of -0.5~ -1.0. In 
fact, in our calculations, we take 5 = —0.5. Since we are mainly discussing the 
significance of FSI, the choice of S does not influence our qualitative conclusion at 
all. 

(ii) With the final state interaction. 

Here, as discussed in our previous work on the FSI, one only needs to calculate 
the absorptive part of the diagram. According to the Cutkosky rule, for getting the 
absorptive part of the loop shown in Fig.l, there are two ways to make cuts, i.e. (1) 
let the V and P be on mass-shell while retaining the t-channel pion off-shell; (2) let 



the pion and P be on shell while leaving V off-shell. The second way refers to a three- 
body decay process and numerical computation for similar triangle diagrams in pp 
annihilation [T5| shows that it is much smaller than (1), so we omit the possibility 
in our later formulation. 

(a) For D + -> K°*ir + -> K°p+. 

In the CM frame of D + where v = (1, 0), the calculation can be greatly simplified. 
To obtain the absorptive part of the loop, for example T3, one can just start from 
eq.(l), replace e&, by ±0) 2 <KPi " ™D 5 (pI ~ ™l)Hto" + P -^-) and add the 
effective vertices of the strong interaction as well as the propagator of the t-channel 
pion. 



x {-2( Cl + ±-(1 + 5)c 2 )ujM~ D [{^Ez - \{U - Z 6 )(M* -m{- ml))(-g , + 
(6 + i^i-g^ + WgL) + i^ k Vimk {-g ltl + 

+ [2(c 2 + ^(1 + 6) Cl )f K , mi 4M~ D {U-%, + ^) - + ^r^))]}- 

(5) 

(b) For L>+ -> ^°p + -> /i"*°vr + . 

The amplitude T4 has a similar form as T3 with some changes, and can be easy 
to obtain from |3|. Saving space, we do not give the explicit expression here. 

In the expressions, for simplicity of bookkeeping, we set pi = p K „ , p<i = p K , p% = 
p n +, p 4 = p p and mi = m K „ , m 2 = m K , = m n +, m 4 = m p . 

F{p 2 K o) is an off-shell form factor for the vertices pirn and K* K°tc. Because we 
use the experimental data where all the three particles are real and on mass-shell 
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to fix the effective coupling at the vertices, in the case of the pion being off-shell, a 
compensation form factor is needed and it is 



A 2 - ml 



\2 



(6) 



with A in the range of 1.2 ~ 2.0 GeV. 

From the above equations and employing the helicity-coupling amplitude formal- 
ism given by Chung [IB|, the whole calculation is straightforward though tedious, 
and a direct comparison at the amplitude level is feasible. We present our numerical 
results in the next section. 



III. The numerical results. 

The two strong coupling constants g pn7t and g K * Kn can be obtained from the p and 
K* decay width, respectively. From the newest PDG data [17], we have g p7T7T = 6.1 
and g K , Kn = 5.8. The values of c\ and C2 are taken from Ref . [O , i.e., c\ = 1.26 and 
c 2 = -0.51. 

As aforementioned, we take the non-factorization factor S to be -0.5. The £j's 
(i=l,...,6) have simple Gaussian forms or polynomials. Their explicit forms and 
parameters can be found in ref. j|] . There are three sets of parameters for aj and b{ 
(notation in ref. 0), which look very different. We substitute all the three sets into 
our expressions to carry out the calculations and compare their results. 

(i) First, we calculate the decay rate without taking into account the FSI, i.e. 
we only use eqs. (0) and (0). The results are listed in Table I. 





Fit 1 


Fit 2 


Fit 3 


Exp. 


T{D+ -> K°p + ) (xlO- 13 GeV) 


1.0 


1.0 


0.65 


0.41 ±0.15 


T{D + -> K*°7i + ) (xlQ- rs GeV) 


0.09 


0.049 


0.079 


0.124 ±0.025 
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Table I. Results without considering FSI and with three sets of parameters from M . 



Here Fitl, Fit2 and Fit3 correspond to three different sets of parameters for the 
functions & (i=l,...,6) from Roberts et al.[§ fitting D -> K^lu; D -> K^lu, B -> 
K^J/ip and D — > K^Hu, B — > KJ/tp, B — > K^)', respectively. One notices that 
even though the parameters in the three fits are quite different, the results are rather 
close to each other. These results are still qualitatively consistent with that obtained 
by authors of ref.|J without the FSI, namely the calculated rate for D + — > K*°tt + 
is lower than the experimental value by 1.4~1.6 times, while for D + — > K°p + , the 
calculated number is 1.6 ~ 2.5 times larger than data. 

(ii) The FSI contribution. 

Since we only consider the absorptive part of the loop, we can get the imaginary 
part of the FSI amplitude. The dispersive part can be obtained by the dispersive 
relation with some cut-off parameters and should be of the same order as the ab- 



sorptive one fll5| . The absorptive part of the amplitude gives a lower bound on the 
effect of FSI. 

For a clean comparison and to avoid ambiguity, we use helicity amplitudes. In 
the helicity picture, for the CM frame of D + all momenta of the outgoing P and V 
are aligned along the z axis (or oppositely). Even though the mesons in the loop are 
real, their momenta can deviate from the z axis by an angle 9. So in this scenario, 
for a reaction D — > PV, only e*(0) polarization of the final outgoing vector meson 
contributes. In Table II, we present the numerical value for the absorptive part of 
the FSI loop with A = 1.6GeV, and as a comparison also the helicity amplitude 
without FSI evaluated in the same theory. 
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Fitl 


Fit2 


Fit3 


T(D+ -> K"p + -> K* u 7r+) 


-i0.220 


-i0.223 


-i0.178 




-i0.071 


-i0.051 


-i0.062 


T(£>+ -> 


0.52 


0.52 


0.42 


T(£>+ -> K*°n + ) 


0.15 


0.11 


0.14 



Table II. 



In the table, we only keep the relative values of the helicity amplitudes with and 
without the FSI, while dropping out any common factor such as (G f / V^)V cs V* d etc. 
The values given in Table II are calculated with A = 1.6 GeV. In fact, our numerical 
computation show that as A varies from 1.2 to 2.0 GeV, the corresponding results 
in the first two rows of Table II can change by a factor 2. As one takes a more 
restricted region for A, the results are not very sensitive to the A value. 

Obviously, the FSI effect of T(D+ -> K°p+ -> J?*°7r+) is stronger than that of 
T(D + — ► K*°7i + — > K°p + ); we will discuss the results below. 

IV. Conclusion and discussion 

Above, we calculate the decay width of D + —>■ K*°tt + and D + —>■ K°p + in terms 
of the HQET and consider the contribution of the t-channel single pion exchange to 
the inelastic scattering K*°tt + <-> K°p + . 

From Table I, we notice that the directly calculated value without considering 
theFSI for D + — > K*°-k + is lower than the experimental value by about 1.5 times 
while for D + — > K°p + it is 1.6~2.5 times larger. In our calculations, three different 
sets of parameters from Roberts et al.[|] are used for the functions & (i=l,...,6) in 
the Gaussian forms (or polynomials). The three sets of parameters are obtained by 
fitting (1) D -> KWlu; (2) D -> K^lu, B -> K^J/ip; and (3) D -> K^lu, B -> 
KJ/tp, B — > Kip', respectively. Obviously, the third set of parameters fits more sets 
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of data, and happens to give the best fit (Fit3) to the data in our cases. But there 
is still a discrepancy. 

Table II shows that the FSI effects are of similar order of magnitude to the 
direct production rates. Especially, the amplitude T(D + — ► K°p + — > K*°7T + ) is 
quite large. This will reduce the rate for the K°p + final state and raise the rate 
for the K*°7T + final state; therefore it will bring the theoretical results into better 
agreement with the experimental data. 

In this work we are not trying to fit the data numerically. We do not introduce 
any free parameter. We just check what is the theoretical prediction for the processes 
D + — ► K*°7i + and D + — > K°p + , based on the more reliable and widely accepted 
theoretical framework |4]] without considering the FSI. We find that the calculated 
values without considering the FSI obviously deviate from the experimental data. 
Then we investigate the magnitude of the FSI effects due to single pion exchange 
inelastic scattering with the unitarity approximation which has been testified in 



many practical processes and proved to be reasonable || [TI|. Our results indicate 
that the single pion exchange process K°p + <-> K*°tt + has significant effects for the 
processes D + — > K°p + and D + — > K *°7r + , and may be the reason for the discrepancy 
between the experimental data and earlier theoretical predictions. 

We cannot give more accurate results, since the dispersive part of the inelastic 
scattering amplitude is very model-dependent, so in this work we are not going to 
adjust the cut-off parameters for fitting data to cause some mess and uncertainty at 
this stage. 

Our discovery of the significance of the single pion exchange inelastic scattering 
in the FSI may have important applications to many other D and B decays. For 
example, for B — > DK, Zheng J7| concluded that the FSI from elastic scattering is 
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negligible. It is noted that for the elastic scattering the lightest exchanged mesons 
are a or p, while the inelastic scattering D*K* <-> DK can be realized by exchanging 
pions and may give more significant contributions. The single pion exchange inelas- 
tic FSI may play an important role in the channels which are relevant to evaluating 
CP violation in some channels and precise measurement of the CKM matrix entries. 
The investigation of those effects is in progress. 
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Figure 1: The single pion exchange inelastic FSI loop for 
D+ -> p+K° 



K*°ir+ and 
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